Josephson vortex interaction mediated by cavity modes: 
Tunable coupling for superconducting qubits 
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A quantum-mechanical model for the interaction of Josephson vortices (fluxons) embedded in su- 
perconducting transmission line is presented. The vortices interact through emission and absorption 
of linear waves (electromagnetic cavity modes). We show that in a classical regime this peculiar 
type of interaction is determined by the product of instantaneous velocities of fluxons. In a quantum 
regime, this property provides tunable coupling between vortices which can be used for entanglement 
of vortex qubits. The physical mechanism of the vortex interaction is similar to that proposed for 
qubits based on trapped ions. Different types of transmission lines mediating the vortex interaction 
are proposed. 
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PACS numbers: 03.67.Lx, 03.75.Lm, 74.50.+r 

It is well established that Josephson vortices (magnetic 
fluxons) can be trapped in spatially extended Josephson 
structures 0, 0- Moreover, in the presence of an ex- 
ternally applied magnetic field and dc bias these pecu- 
liar nonlinear objects propagate along a system or are 
trapped around a particular point 0, 0, S ■ I n a clas- 
sical regime the vortex motion is mapped to the classical 
mechanics of a macroscopic particle. As we turn to a 
quantum regime realized at low temperatures and for a 
small energy of the vortex in an external potential, the 
vortex is predicted to behave quantum-mechanicahVe. 
g. show tunneling effect and discrete energy levels [l,Uj]. 
Quantum tunneling of vortices has been observed in dis- 
crete arrays of small Josephson junctions [8|, and recently 
quantum tunneling of a single vortex was observed in 
long annular Josephson junctions 9]. A quantum state 
of a macroscopic vortex in an artificially created poten- 
tial can be used in order to implement a particular typ e 
of superconducting qubits named vortex qubits |lCt . 

Up to now, two types of vortex qubit prototypes have 
been investigated. The first prototype is the heart- 
shaped Josephson junction shown in Fig. ^i, for which 
rather simple state preparation and reado ut p rocedures 
have been already verified experimentally pT| . In these 
qubits the external magnetic field is uniform and the 
spatially-dependent potential profile for a vortex is tai- 
lored by shaping the long junction into a heart form. The 
second vortex qubit idea is based on creating a desired 
potential for the vortex by using local magnetic fields 
|12| generated by control currents I\ and I2, see Fig. 
Both cases can be described by a local double-well po- 
tential U(yi) for the vortex, see Fig. QJ:. Here, y^ is the 
coordinate of the center of the vortex. 

The qubit shown in Fig. QJj has the advantage that 
its state preparation and readout can be done by us- 
ing Rapid Single Flux Quantum (RSFQ) circuits. RSFQ 
logic is the natural choice for an interface between vor- 
tex qubits and room-temperature electronics. The field- 
injector based vortex qubits with RSFQ preparation and 
read-out interfaces have been also recently tested in the 
classical regime ^3|. The entanglement that persists 
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FIG. 1: Possible realizations of vortex qubits. a) Heart- 
shaped long Josephson junction placed in the magnetic field 
H [ill . b) Long Josephson junction with local magnetic field 
injectors, c) Double-well potential for a vortex in the above 
two systems. Two vortex states are indicated as A and B. 



in the quantum regime in a multi-qubit circuit should 
be possible to read out by RSFQ circuitry. The vortex 
qubits are naturally scalable to larger quantum circuits. 

In a system containing a number of vortices they in- 
teract with each other. The simplest type of the inter- 
action, namely direct static magnetic interaction, occurs 



when the vortices are located in the same long Joseph- 
son junction, see Fig. In this case the interaction is 
due to the presence of a mutual inductance between the 
two vortex locations and is determined by overlapping of 
the "tails" of vortex magnetic fields. Such an interaction 
decays exponentially [lj] with the distance between the 
vortices on the characteristic scale of the Josephson pen- 
etration depth Aj and depends on their coordinates. In 
the quantum regime the coupling between vortices leads 
to a quantum entanglement, which is the necessary con- 
dition in order to realize a quantum computation. 

The static interaction exponentially depends on the 
distance between vortices so that only coupling between 
the nearest neighbors has to be taken into account. In 
coordinate representation the static interaction term be- 
tween two vortices can be written as Uq exp(|yj — yj\/\j), 
where yi, yj are the coordinates of the interacting vor- 
tices, and Uq is the characteristic interaction energy. If 
we switch to the representation of wave functions corre- 
sponding to the two energy levels of single vortices, the 
interaction term in the Hamiltonian can be written as 
H m t = ^2,Uq exp(|aj — a,j\)a\a 3 z , where <2j, aj are the av- 
erage coordinates of corresponding vortices normalized 
to Aj. Note here, that with this type of interaction the 
chain of vortices is mapped to a 1-D Ising model of in- 
teracting spins 1/2 or to a system of nuclear spins |l5j . 

A severe drawback of the static type of interaction is 
that it can not be turned off and addressing of a chosen 
qubit (pinned vortex) will immediately affect all other 
qubits. Thus, in this realization it is easy to obtain en- 
tangled states but it remains a hard problem to address a 
single vortex, e.g. switch from the ground state to the ex- 
cited state (|0 >^ |1 > transition) of a particular vortex 
without affecting others. 

In order to avoid the above problem, we propose for 
coupling of vortex qubits placing them in a supercon- 
ducting transmission line (Fig. [2h and 12k). In this case 
the vortices display a peculiar indirect dynamic inter- 
action through virtual excitation (absorption) of linear 
electromagnetic waves (EW) propagating in a transmis- 
sion line. This interaction can be realized at least in 
two different configurations. One is a lateral configura- 
tion shown in Fig. |2Jd, where two or more long Joseph- 
son junctions are incorporated in a common system with 
a superconducting ground plane. In such a system the 
superconducting electrodes and the ground plane pro- 
vide a low dissipative transmission line for propagation 
of EW. Thus, a moving (oscillating) vortex excites the 
EW (emission of photons) which propagate through the 
transmission line and interact (by absorption of photons) 
with other vortices. Note, that this principle allows to 
realize a quantum-mechanical interaction similar to the 
one proposed for chains of trapped ions [16l Il7| . In the 
latter system single-qubit operations address ion energy 
levels but the interaction is just due to the excitation of 
phonons in a chain of trapped ions. Similar realization 
of couplin g vi a lumped L-C circuit has been proposed 
earlier |l8lll9| for superconducting qubits. 



vortex qubits 1 



(a) 




control line 



qubits 1 2 3 



(b) 




(c) 







ground plane 




passive stripline 


EW 




" — s* ■ t 


77~ 





long junction 



FIG. 2: Sketch of long Josephson junction circuits allowing 
for different types of vortex interaction, a) Static coupling 
between vortex qubits can be realized by placing them side 
by side in a long junction. The distance d between qubits 
has to be larger (but not much larger) than Aj. b) The 
dynamic coupling between vortex qubits located in different 
junctions. The junctions are placed above a superconducting 
ground plane, c) The dynamic coupling between vortex qubits 
separated by large distance in a single long junction. Note a 
passive stripline behind the junction. It provides the coupling 
between qubits via electromagnetic waves (EW) propagating 
in the stripline. 



Anther possible realization of the dynamic coupling is 
a configuration with far-separated vortices located in the 
same Josephson junction. These vortex qubits indirectly 
interact with each other through the EW excitations in a 
passive transmission line inductively coupled to the junc- 
tion, see Fig. [2;. Similarly to a lateral configuration, the 
EW propagating in transmission line provides a fast and 
low dissipative source of interaction. 

In the following, we quantitatively analyze the lateral 
configuration of vortex qubits located in different long 
junctions placed above a superconducting ground plane 
(Fig. Ob)- A classical dynamics of such a system is deter- 
mined by time-dependent Josephson phases ip(y — yi,t), 
where y is the coordinate along Josephson junctions and 
yi is the coordinate of the center of the vortex located in 
the i-th Josephson junction. The EW propagation along 
the transverse coordinate x is described by a set of equa- 
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where the coordinates of Joscphson junctions, a 

is the width of the Josephson junction in x- direction, 
V(x, t) is the voltage between the superconducting elec- 
trodes and the ground plane, I{x, t) is the current flow- 
ing in x direction along the superconducting electrodes. 
The parameters L (the inductance of superconducting 
electrodes per unit length) and C (the capacitance be- 
tween superconducting electrodes and the ground plane 
per unit length) determine the properties of a transmis- 
sion line created between the superconducting electrodes 
and the ground plane. The parameter R s describes the 
dissipation due to the surface losses (in general, they are 
frequency dependent) in the transmission line. For not 
very high frequencies these losses in the superconducting 
line are very low. Notice here, that without the term 
depending on the Josephson phases this set of equation 
describes the propagation of a weakly decaying linear EW 
with the spectrum ui(k) = ck, where c = 1 / \/ LC pcj . 
In the presence of vortices we obtain the wave equation 
where the vortices are the local sources of perturbation: 
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Thus, we obtain a time Fourier transformation of the ac 
current distribution in the transmission line I u (x) in the 
form: 

Ux) = C—J2oj 2 G u (x,Xi) [ dte iut yi (t) , (3) 
ae ^ J 



where G u (x, Xi) = G^ix — Xi) is the Green function of 
the linear equation: 



d 2 G u} {x 1 x l ) 
dx 2 



The equation of motion for a vortex center of mass in 
the presence of a local double- well potential U(yi) illus- 
trated in Fig. ^ is written as 



jji + 



2e\jdU{ Vl ) I( Xi ,t) 



hi, 



<:■() 



dyi 



l cO 



= 



(5) 



where u p and I c q are the plasma frequency and the crit- 
ical current of a Josephson junction (in the absence of 
vortices), respectively 0, 0- Thus, the interaction en- 
ergy between vortices can be found as: 
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By making use of 10 and (@J we obtain 

v)' 
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X E/ dt m(t)yj(h) J dueW-^Gufa-xj) . (7) 

This expression can be symmetrized over the indexes i 
and j as 
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where the kernel of interaction 

K(t,x) -- i / duuje-'^GuW.,] 



(9) 



depends on the parameters L and C in our model. In the 
limit of a small propagation time (r = \xi — Xj\V LC <C 
T osc , where T osc period of slow vortex oscillations between 
two quantum states) we neglect a delay in the wave prop- 
agation, and obtain the classical interacting potential be- 
tween vortices located in the Josephson junctions i and 
3- 
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(10) 



where the constant of interaction Aij is determined by 
the properties of the transmission line given by For 
a simplest case of a transmission line described by © we 
obtain 
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Thus, in the case of an indirect dynamic coupling the 
vortex interaction is determined by a product of vortex 
velocities. 

In the quantum regime this dynamic interaction is de- 
termined by operators of vortex momentum (pi and pj): 
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where m/ is the effective mass of a vortex |lOjj . The total 
Hamiltonian H of indirectly interacting vortices can be 
written as 

If a single vortex potential U(yi) has a double-well 
form, the quantum-mechanical description of a single vor- 
tex can be reduced to lowest energy levels in each well. 
These levels are characterized by the energy difference 
~ B z and the tunneling amplitude ~ B x between states 
[To) . The parameters B z and B x are determined by the 
properties of a single vortex potential U(yi). Because the 
average value of the momentum operator < pi > is zero 
in these states, we obtain that the interaction potential 
is determined by tunneling amplitudes B x (i). By mak- 
ing use of the spin representation of the Hamiltonian we 
write 

H = ^[fl,(i)^+B a (i)^]+i^fl,(i)B s ^ii 3 -5-*^ 

i i^j 

(13) 



A great advantage of this dynamic type of interaction 
in respect to a direct static one is that the former one 
is absent when both vortices are in the ground state. 
With this type of Hamiltonian it is possible to address 
a single vortex or a bunch of vortices, if that is required 
|l9j . Indeed, by turning off the tunneling amplitude B x 
of a particular vortex we exclude it from the computation 
process, but other vortices can still be manipulated. 

In conclusion, we proposed here and quantitatively an- 
alyzed an interaction of vortices through an exchange of 
EW propagating in a superconducting transmission line. 
In a classical regime, this interaction depends on the ve- 
locities of vortices. As we turn to a quantum regime, this 
type of interacting allows for a tunable coupling between 
vortices that is a crucial condition for implementation of 
such vortices as superconducting qubits. The addressing 
of these qubits can be provided by a particular combina- 
tion of dc pulses allowing to change the parameters B z 
and B x . 
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